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Abstract 



We analyze the Einstein equations for a spherically symmetric static distri- 
bution of matter which satisfies a given constitutive relation for elastomechanics. 
After reducing the equations into a system of Fuchsian ODE for certain scalar 
invariants of the strain, we show that, for a given constitutive relation and a given 
value of central pressure satisfying certain compatibility conditions, there exists a 
unique regular solution near the center. In the case when the constitutive relation 
is given by a quadratic form of strain, we further show that the solutions stay 
regular up to the boundary of the material ball. 



f Current Address: 134 Ivy Dr. Apt.5, Charlottesville, VA 22903, U.S.A. 



§1 Introduction 



It is common in relativistic astrophysics to regard a star as a massive ball con- 
sisting of perfect fluid. Assuming that the ball is static, the corresponding Einstein 
equations reduce to the well known Tolman-Oppenheimer-Volkoff equation. Vari- 
ous authors have discussed properties of solutions of this equation. The existence 
theorem for this equation under general equation of state has been established by 
Rendall and Schmidt [1]. In the same article, as well as in a more recent article 
by Makino [2] , the issue of fmiteness of the ball and the regularity of solutions up 
to the boundary of the ball has been discussed. 

Anisotropic material balls also have been objects of study and many particular 
solutions of the corresponding Einstein equations have been derived (for example 
[3], [4], [5] and [6]) under various assumptions imposed on the stress-energy tensor. 
But, unlike in the case of the perfect fluid ball, there hasn't been any general 
statement concerning the existence of solutions for anisotropic balls. The purpose 
of this paper is to resolve this problem for the elastic solids. To achieve this 
purpose, we use the relativistic elasticity theory formulated years ago by Kijowski 
and Magli [7]. Although there have been other complete theories of relativistic 
elastomechanics (for example [8] and [9]), the formalism of Kijowski and Magli 
seems to suit better for our purpose. 

The main idea common in the modern theories of relativistic elasticity is that 
a material property is chosen by choosing a constitutive relation — the relation 
between the energy stored at a particle and the strain experienced by that particle. 
Once a constitutive relation has been chosen, the stress can be written as a function 
of strain and hence the material variables are replaced by the dynamic variables 
which describe the strain. Under this scheme, the Einstein equations for the 
spherically symmetric static elastic body can be written as a closed system of 
ODE for a variable describing the material configuration and two variables for the 
spacetime metric. In [7] some families of solutions of these equations have been 
produced, but no theorem regarding the existence of solutions of this system under 
general constitutive relations has been given. 

In this paper, we show that for a given constitutive relation and a prescribed 
value of central pressure which satisfy certain compatibility conditions, there ex- 
ists a unique smooth solution on a neighborhood of the center to the equations 
mentioned above. To prove the local existence and uniqueness, we use the theorem 
by Rendall and Schmidt [1] which states the existence, uniqueness and regularity 
of solutions for a certain class of Fuchsian ODE systems. 

It is also shown that, in the case the given constitutive relation is a quadratic 
form of strain, the solutions can be extended to the boundary (where the radial 
stress vanishes), so that a Schwarzschild vacuum can be continuously joined outside 
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the ball. This is done by showing that the solutions are bounded on any finite open 
interval on which the radial stress stays positive. The most crucial step is to show 
that the mass-radius ratio is bounded from above by a constant < -|. The similar 
issue actually arises in many other matter models under spherically symmetric 
setting. For perfect fluid and anisotropic fluid, this estimate has been made by 
Baumgarte and Rendall in [10] under the assumption that the tangential stress (in 
the case of anisotropic fluid) and the equation of state relating the energy density 
and the radial stress have been given and are sufficiently regular. In our paper the 
same estimate is made but in a different context; the regularity of tangential stress 
is not assumed a priori, but it is a consequence of the mass-radius ratio estimate. 
For Vlasov field, with isotropic or anisotropic stress, the same estimate has been 
established by G. Rein in [11]. We establish the bound of mass-radius ratio for 
the elastic material by adopting the method used by Rein in [11]. 

Throughout the paper, we use the spacetime metric signature (—,+,+,+) 
and the units c = G = 1. 
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§2 A Brief Review of Relativistic Elasticity and 

Generalized Tolman-Oppenheimer-Volkoff Equations 

In relativistic continuum mechanics, it is common to describe the configura- 
tion of particles by a ma,p(material map) from a spacetime into a three dimensional 
Riemannian manifold( material space). This map £ : (M ,gij) — > (X s , 7a&) must 
satisfy the following conditions: (i) its differential D£ has full rank everywhere and 
(ii) at each p e M, the solution w l e T p (M) of the equation u 1 ^ = 0,6 = 1,2,3 
is timelike. Both the projection rjij := gij + UiUj of the spacetime metric along 
u l and the pull-back (£*7)y := i^jlab of the material metric act as Riemannian 
metrics on the subspace of T P M which is orthogonal to u l . When X represents an 
elastic body, the material metric 7^ describes the infinitesimal distance between 
particles in a locally relaxed state. So, we say that an event p G M is at a re- 
laxed statef if these two tensors agree on the orthogonal subspace in T p M of u l . 
There are many different ways of describing the deviation of 77^ from (£*7)y and 
any of them can be called a relativistic strain tensor. Among them we choose the 
following definition of strain as in [7] . 



Note that the operator to which "log" is applied is positive and self-adjoint, so the 
strain is well-defined by (1) for any choice of 7 a & and £. One can also show that 
this tensor is symmetric and is spatial in the sense that SijU^ = 0. 

The main idea in the constitutive theory of relativistic elasticity is that the 
strain tensor should fully determine the elastic energy E stored at each particle, 
that is, E should be a function of Sij. We call this function a constitutive relation 
for the material. Choosing this function E = E(Sij) is equivalent to choosing a 
particular type of material to study. 

The elastic energy E as a function of strain determines the stress-energy tensor 
which satisfies the conservation law E>{I l - = 0. The relationship between the 
stress-energy tensor and the strain tensor generalizes the stress-strain relationship; 
the details of its derivation can be found in [7]. We will quote the result for an 
isotropic elastic material only. 



f But, for a body under extreme pressure, a relaxed state might not exist. See 
[8] for a detailed discussion and an alternative formulation. 




(1) 
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where "tilde" means "traceless part of" , n is the particle number density measured 
with respect to the pull-back metric an d 

a := trS (3) 
:= ±tr[S • S] (4) 

9:=±tr[S-S-S] (5) 

The quantity a measures the compression rate of the material relative to the 
relaxed state. More precisely, the particle number density n with respect to the 
projected spatial metric rjij is given by 

n = n e~ a . 

The energy density is therefore p = nE = n e~ a E. The part between the round 
brackets in (2) (trace-free part) vanishes identically if E depends only on the 
compression a, i.e. when the material is a perfect fluid. 



Now we will consider the case when the spacetime metric gij is spherically 
symmetric and static. We will assume the following form of spacetime metric. 

ds 2 = -e v dt 2 + e x dr 2 + r 2 d9 2 + r 2 sin 2 6># 2 

Here the functions v and A are supposed to be smooth on M and depend only on 
the radial coordinate r. This is not the most general form of spherically symmetric 
spacetime metric since in general e~ x can vanish for some r > 0, but we will see 
in §4 that this cannot happen in the interior of material ball. We choose X = R 3 
as the material space and the Euclidian metric f 

dy 2 + y 2 d6 2 + y 2 sm 2 9 # 2 

as the material metric. We also assume that the particle distribution at the relaxed 
state is homogeneous, i.e. n a is a constant. Under this assumption, we may set 
n Q = 1 without loss of generality. A material map is given by specifying three 
functions y(t, r, 9, 0), 9(t, r, 9, (p) and r, 9, (p). But, by the assumption that the 
spacetime is static and spherically symmetric, y must depend on r only. After 
adjusting the angular coordinates for X if necessary, we may further assume that 

f For more general material metrics, see Appendix. 
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9 = 6 and = 0. So, a material map is determined by specifying only one function 
y : [0, oo) — > [0, oo). We will allow only smooth material maps, so y(r) = rz(r 2 ) for 
some smooth function z : R — > R. We also insist that y'(r) > for all r G [0, oo) 
in order to ensure that the material map has full rank everywhere. 

With a function y(r) given for the material map, we can now form the strain 
tensor from (1). We have 

S = -log 4=== ( e^dr 2 ) - log^ (r 2 d9 2 ) - log^ (r 2 sin W) (6) 

v / e A(r) V / r r ' 

To write the stress-strain relationship (2), we need the partial derivatives of the 
energy with respect to the invariants of the strain. As a consequence of spherical 
symmetry, there are only two independent invariants of the strain tensor. Hence, 
the elastic energy E can be considered as a function of these two independent 
invariants of the strain. It turns out that the following choices of two invariants 
make the expression of the stress-energy tensor much simpler. 

v ■= a + u = -31og (J^j (8) 

The invariants a, 8, and 6 are related to u and v as follows. 

1 2 

a = v-u, 8 = -u 2 6 = — -u 3 (9) 

3 27 K ' 

Using (2) and (9), we now get the stress-energy tensor. 



rp. . _ p U~V 
± lJ — V 



Ee v dt z + ^-e x dr l +[-£-- -(^ + r\d9 2 + s\n 2 6d<p 2 ) 

au \ au 2 ou ov ' 



(10) 



Let us pick notations for the components of the stress-energy tensor. 

p := e u ~ v E (12) 

P:=e u-v^± ( 13 ) 



e u ~ v E 




e u- v dE 




du 




2 [ du 


+ ^» 




3 dE 



dE 



(14) 
(15) 
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The Einstein equations can be now written. 



8 7 rr 2 p = e" A (rA / -l) + l (16) 
8wr 2 P = e~ x (ru' + 1) - 1 (17) 



Svr Q = \e~ X (v" + + V - A') - ±i/A') 



(18) 



It should be remarked at this point that this is a closed system for three unknowns 
{y, A, v} since E is a function of u and v, which are again functions of y, y', and 



Next, we will modify the Einstein equations (16)-(18) into the generalized 
Tolman-Oppenheimer-Volkoff equations by introducing a new variable 

w := ^ (1 - e" A ) (19) 

With this new variable, the equation (16) becomes 

rw' = 4yrp - 3w (20) 

The other two equations become 

P' = - r(l - 2r 2 w)~ 1 (4nP + w) (P + p) (21) 

r 

and 

v' = 2r(l - 2r 2 w)- 1 (47rP + w) (22) 

The equation (22) can be integrated after all the other solutions have been found. 
So, we will focus on the equations (20) and (21) only. 

We will assume throughout the paper that the constitutive function E(u, v) 
is smooth. We will be also interested in only smooth solutions, that is, smooth 
material map and smooth metric coefficients. With these assumptions, it follows 
from (12)-(15) that all the variables which appear in the equations have to be 
smooth. So, we may consider them as smooth functions of x = r 2 defined on R. 
The equations (20)-(21) then can be expressed using the new independent variable 
x = r 2 as follows. 

2xP' = -60-x(l - 2xw)- 1 (4nP + w)(P + p) (23) 
2xw' = Anp - 3w (24) 
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Here, and throughout the rest of the paper, the "prime" indicates the derivation 
with respect to x. We also have an auxiliary equation from (7) and (8), 

2xv' = 3-3(1 -2xw)~h u (25) 

The equations (23)-(25) form a closed system for the variables u, w, and v. 
If a set of smooth solutions {u(x),w(x),v(x)} are found, then the solutions for 
the original equations will be determined algebraically or by a simple integration. 
The function y will be determined from (8) since it can be explicitly solved for y 
as y(r) = re~ v ^ r \ The metric variables A and v will be determined by (19) and 
(22) respectively f 



f (19) can be solved for A only if we know that w < (2r 2 ) -1 . It will be shown 
in §4 that indeed this inequality holds on any interval on which the radial stress 
P never vanishes. 
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§3 Local Existence and Uniqueness of Solutions 



This section is concerned with the integrability of the equations (23)-(25) on 
a neighborhood of x = 0. After suitable substitutions and rearrangements we can 
write the equations explicitly in terms of u, w and v. 

9 , _ 3(£ M - E uv )(l- Vbe u )- 3(E U + E v ) xb(4ne u - v E u +w)(E+E u ) 

(E + E ) E +E [ ' 

\ J -'u ~ J -'uu) J -'u ~ i-'uu 

2xw' = 4ne u ~ v E - 3w (27) 
2xv' = 3 - 3Vbe u (28) 

where 

b:=(l-2xw)~ 1 (29) 

and E is a smooth function of u and v. These equations are singular at x = 
and therefore the existence theorem for regular equations cannot be applied here. 
But, there exists a theorem which states existence, regularity and uniqueness of 
solutions for this type of equations, which we state below. 



Theorem 1 (Rendall and Schmidt). Let V be a finite dimensional vector 
space, N : V — > V a linear map all of whose eigenvalues have positive real parts, 
and G : V x (— e, e) — > V and g : (— e, e) — > V smooth maps, where e > 0. Then, 
there exists 5 < e and a unique bounded C 1 function f : (—5, 0) U (0, 5) — > V which 
satisfies the equations 

x %,+ N -f = xG ( x i /(*)) + 9(x) (30) 

Moreover, f extends to a smooth solution of (30) on (—5, 5). If N, G and g depend 
smoothly on a parameter t and the eigenvalues of N are distinct, then the solution 
depends smoothly on t. 

Our equations are not in the form (30) yet. For instance, the equation (26) 
contains terms which are non-linear in u and v and those terms do not have a factor 
of x. One way to introduce a factor of x in those terms is to write u = u Q + xu\ 
and v = v a + xv\ and to rewrite the equations for the new variables u±, w and v±. 
By the way, (28) implies that v(0) can be chosen freely and u(0) = necessarily. 
So, we let 

U = XUi, V = V + XVi 



9 



Then the equations (26)-(28) become 

3(E U - E uv )(l - Vbe u ) - 3(E U + E v ) 



2xu'± + 2u\ 



x(E u + E uu ) 
b(4ive u - v E u + w)(E + E v 



E u -\- E uu 

2xw' + 3w = 4ne u - v E 

2xv[ + 2v 1 = 3a;- 1 (1 - Vbe u ) 



(31) 

(32) 
(33) 



We will now show that, if E(u, v) and v a satisfy certain conditions, then these 
equations can be rearranged into a system of the form given in Theorem 1. Let 
U := (u!,w,vi). 

Lemma 1. Suppose (i) E u (0, v ) + E uu (0, v Q ) ^ and (ii)E u (0, v) + E v (0, v) = 
for all v. Then, the equations (31)-(33) can be rearranged into the form 

2xU' + AU = xG(U,x) + V (34) 

Here Q : R 3 x R — > R 3 is a smooth map, V is a constant vector in R 3 and 





c 


k 


o\ 


A = 





3 







\3 


3 


2/ 



where k is a constant which depends on vq- 

Proof. Throughout the proof, "O(x)" will represent a generic function of the 
form xf(ui, w, Vi, x) where / is a function which is smooth on a neighborhood of 
the set x = in R 4 . 

(1) The First Equation: The first term on the right hand side of (31) is 
3(E U + Euu)' 1 [{Eu - E uv )x- 1 (l - Vbe u ^j - x~\E u + E v )j 

We have 

Vbe u ) =x~ 1 1 - (1 -2xw)-h XUl j 

= x" 1 [1 - (1 + xw)(l + xux) + 0(x 2 )] 



x 



= —w — ui + 0(x) 
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(35) 



and 

(E u - E uv ){u, v) = (E u - E uv ){0,v o ) + 0{u) = (E u - E uv ){0,v o ) + 0{x) 

So, 

(E u - E uv )x~ 1 (l - Vbe u ) = (E u - E uv )(0,v o ) ■ (-w - Ul ) + 0(x) (36) 

Next, using the second assumption, (E u + E v )(0, v) = 0, we have 

(E u + E v )(u,v) = x Ul (E u + E v ) u (0,v o ) + xv t (E u + E v ) v (0,v o ) + 0(x 2 ) 
= xui(E uu + E uv )(Q,v ) + 0(x 2 ) 

So, 

x- 1 (E u + E v )(u,v) = u 1 ■ {E uu + E uv )(0,v ) + O{x) (37) 
Combining (36) and (37), we get 

(E u - E uv )x~ l (l - Vbe u ^j - x~\E u + E v ) 

= (E u - E uv )(0,v o ) ■ (-Mi -w)- (E uu + E UV ) (0,v o ) ■ ui + 0(x) 
= -(E u + E uu )(0,v o ) ■ u x - (E u - E uv )(0, v )-w + 0(x) 

On the other hand, by the first assumption, (E u + E^) -1 is a smooth at (0, v ). 
So, (Eu + E^)- 1 = (E u (0,v o )+E uu (0,v o ))- 1 + O(x). Therefore, 

3(E U -E UV )(1- Vbe u )-3(E U + E v ) ( (E u - E uv )(0, v ) \ , nf , 

- —6Ui — 3 — — r W + 0(X) 



x(E u + E uu ) " \(E u + E uu )(0,v o ) 

(38) 

The second term on the right hand side of (31) is of the form b ■ g(u,w,v), 
where g is a smooth function which is linear in w. Since b = (1 — 2xw)~ 1 = l+0(x) 
we have 

b ■ g(u,w,v) 

= g(0, 0, v ) + xu!g u (0, 0, v ) + wg w (0, 0, v Q ) + xvig v (0, 0, v Q ) + 0(x) 

= d + C 2 w + 0(x) (39) 

for some constants C\ and Ci which depend on vq. Therefore, adding (38) with 
(39), the equation (31) becomes 

2xu[ + 2m = -3ui -kw + 0{x) + Vi 
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or 



2xu[ + 5ui + kw = 0(x) + Vi 



for some constant V\. 

(2) Second and Third Equations: Since e u ~ v E(u, v) is smooth, we can write 
this as C\ + XU1C2 + XV1C3 + 0(x) = C + 0(x). So, the second equation is already 
in the desired form. The third equation has been already handled by (35). 

[End of the Proof of Lemma 1] 

Remark: The first hypothesis of this lemma restricts the choice of data v a when 
a constitutive relation is given. But, unless E u + E uu = on the line u = 0, there 
will be always an open interval I such that the hypothesis is satisfied for all v Q G /. 
This tells us that for almost all the choices of constitutive relation, there will be 
enough room to pick up vq with which the first hypothesis is satisfied. 



A local existence theorem immediately follows from Theorem 1 and the above 
lemma. 

Theorem 2. Assume E(u,v) and v Q satisfy the hypotheses of Lemma 1. Then 
there exists a positive constant S such that the equations (27)-(28) admit a unique 
set of C°° solutions {u,w,v} on the interval (—5,5) satisfying v(0) = v Q . 

Suppose that the constitutive relation is given in the following form 



where m represents the rest mass energy of a particle and A, B are nonnegative 
constants. Then, the condition (i) is satisfied for all v Q except for v a = 1 + |^ and 
the condition (ii) is also satisfied. 




— (u — v) 2 H 

2 v ' 3 




(40) 



12 



§4 Global Properties of Solutions for Quasi-Linear Models 



In this section we study the global properties of solutions under the assump- 
tion that the constitutive relation is given by (40). The elasticity theory under 
such a form of constitutive relation is called the quasi-linear elasticity because 
under this constitutive relation the stress-strain relationship becomes linear. Un- 
der the quasi-linear constitutive relation of the form (40), the stress-energy tensor 
components (12)-(15) becomes 



P 

P :-- 
O := 
Q:- 



e 

B t 
— ue 
3 



m+—(u-v) + — u 



^ 2B 

A(u — v) H u 

3 



A(u — v) — — u 



(41) 
(42) 
(43) 
(44) 



One of the properties of solutions in which we are interested is whether the 
radial stress P vanishes at some finite radial distance or not. If P vanishes at 
some finite radial distance r = R, we can consider the sphere r = R as the bound- 
ary of the material body. If, further, the spacetime metric stays bounded inside 
this sphere, then a vacuum Schwarzschild solution may be attached outside the 
sphere to obtain a complete C 1 spacetime. Therefore another important question 
is whether the solution stays finite up to the radial distance where the radial stress 
vanishes. 

This section focuses on answering the last question, but the first question will 
be answered in a special case when the anisotropy constant B is sufficiently close 
to 0. From now on we assume that v(0) = v D < is given. Since it(0) = 0, this 
means that we are assuming that the radial stress P at x = is positive. We will 
show below that if P stays positive on a finite interval then all the variables u, v 
and w stay bounded and the equations (26)- (28) stay regular (except at x = 0) on 
the interval. The standard extension theorem of ODE then would imply that the 
interval of existence can be enlarged, i.e. the interval is not the maximal interval 
of existence yet. 



Bounding w. 

Let [0,R), R < oo, be an interval of existence of regular solutions of our 
equations (23)-(25) and assume P > on the interval. Suppose w were not 
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bounded from above on the interval. Then by continuity of w on the interval, 
it would follow that 1 — 2xiw(xi) = for some X\ G [0,R). Then from the 
equation (23), O would have to have a singularity at x = x\. This contradicts the 
assumption that the solution is regular on [0,R). Therefore w must be bounded 
from above. To obtain the lower bound, we use the equation (20) from which it 
follows that w(0) = 47rp(0) > and 

-^(r 3 w(r)) = A%pr 2 (45) 
dr 

for all r. So, r 3 w(r) > for all r. Therefore, w > 0. Moreover, because -^(r 3 w) 
is strictly positive (since p > 0), w must be bounded from below by a positive 
constant if we restrict the interval to [x Q , R) for some x a > 0. 



Bounding e . 

First we will check whether the integral of p is bounded on finite intervals. 
From (45), we have 4Tvp(r) = r~ 2 -^(r 3 w(r)). Integrating this, we get 

pr pr 

An / p(s)ds = rw(r) + 2 w(s)ds (46) 
Jo Jo 

So, it follows that the integral of p must be bounded on any finite interval since 
we know that w is bounded. Now, we are ready to show that e A is bounded. The 
proof of following lemma is based on an argument that has been used by G. Rein 
in [9]. 

Lemma 2. Let [0, R) be a finite interval on which the regular solution exists and 
P > 0. Then, e x is bounded on the interval. 

Proof. Using the original form of the Einstein equations (16) and (17) and the 
generalized Tolman-Oppenheimer-Volkoff equation (21) we can derive the following 
equation. 

(e > (4«P + m) ) =4«e » (J- - - - _ j (47) 
On the other hand, from (41)-(44), it follows that for a choice of constant 



k > max 
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we have 

\P\ < kp and \Q\ < kp 

and hence 

|0| < 2kp 

So, noting that w > 0, we get the following inequality from the equation (47). 

(e^ (4nP + «;))'< Kp (e^) 

for some positive constant K. By restricting the interval to [xo, R) for some x Q > 
if necessary, we may assume that w is bounded from below by a positive constant. 
Then, the above inequality can be replaced by 

(^^(4ttP + wj)' < Kp (e^(4nP + wfj 
Applying Gronwall type argument to this inequality, we get 

(e^ (4nP + w^j (x) < Cexp ^K J p(s)ds^j 

for some positive constant C. So, using (46) and the earlier observation that w is 
bounded, we get 

e^(4nP + w) < C 

on [0,.R). From here, using again the fact that w is bounded from below by a 
positive constant we get e 2 ^ < C. Now, the desired bound < C follows since 
by the equation (22) v is monotonically increasing. 
[End of Proof] 



Bounding the Rest. 

As before, we assume that P > on the interval [0, R). 

Step I. Showing that v, P and O are bounded from above. It follows immediately 
from the equation (25) that v is bounded from above. We get from (23) 

2(xP)' = 2xP' - 2P < -60 - 2P < 6|0| + 2\P\ <k f p(s)ds 

Jo 
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where, in the last step, we used the fact that \P\ and |0| are bounded by kp for 
some constant k. This gives an upper bound of P since we know that the integral 
of p is bounded. The upper bound of O then follows from 

20 = P - A(u - v)e u ~ v < P + 1 

Step II. Showing that u is bounded from above. If B > 0, we have 

ue u = -|oe" < C 

So, u is bounded from above. If B = 0, then we know that A > 0, so we have 

(u-v)e u ~ v = ^r<C 

Therefore u — v is bounded from above. Since v is already bounded from above, 
it follows that u is bounded from above. 

Step III. Showing that u and v are bounded from below and the equation (26) 
remains regular within the interval [0, R). A lower bound of v comes from the 
equation (28) since we know that \fb := is bounded and that u is bounded 
from above, u is bounded from below since 

e v ~ u P =^A+ u-Av>0 

and v is bounded from below. The singularity of the equation (26) outside x = 
can occur only when E u + E uu vanishes. But, from (40), 

E u + E uu = A{u -v) + ^-u + A + ^- = e v ~ u P + A + ^ > A + ^- > 
So, there is no singularity away from x = 0. 



We have shown so far 

Theorem 3. Let E(u,v) be given by (40). Given a constant vo < 0, let [0,R) be 
the maximal interval on which the regular solution o/(26)-(28) satisfying the initial 
condition v(0) = v Q exists. Then P must vanish in the interval unless R = oo. 

Remark: Note that this theorem applies only to the quasi-linear constitutive 
relation while the local result (Theorem 2) could be obtained for much larger class 
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of constitutive relations. Many of the arguments used to prove Theorem 3 rely 
on special features of the quasi-linear constitutive relation given by (40). But, 
the most crucial step — the boundedness of e A - - could have been established 
for a much larger family of constitutive relations, namely those which satisfy the 
dominant energy condition. Note that under the dominant energy condition, the 
quantity \Q\ is bounded by 2p, which was the most important step in the proof of 
Lemma 2. 



Finiteness of Radius 

Suppose B = 0. Then O = 0, so the material is an isotropic fluid. We will 
show that the radius is finite. Suppose P > for all x. Then, by the monotonicity 
of P (see equation (23)), we know that P tends to a limit as x — > oo. So, P' tends 
to as x — > oo. Then, again by the equation (23), it follows that p — > and P — > 
as x — > oo. But 



We can see from here that it is impossible for p and P both tend to while P > 0. 
Therefore P must vanish somewhere. By the continuous dependence of solutions 
on parameters (Theorem 1), it immediately follows that for all sufficiently small 
B 1 the solutions must represent balls with finite radius. 

§5 Discussion 

We have shown that, given a smooth constitutive relation and a prescribed 
value of the stress at the center, there exists a unique static relativistic elastic 
ball. We also have shown, if the elastic energy is quadratic in the strain tensor, 
then either the radial stress vanishes at a finite radial distance before the solution 
becomes singular or the material fills up the whole space. If the radial stress van- 
ishes at a finite radial distance r = R, then we can join a Schwarzschild vacuum 
outside the sphere r = R and obtain a complete asymptotically flat C 1 solution of 
the Einstein equations. In the case that the radial stress never vanishes, which is 
certainly possible in theory, we can ask whether the solution represents an asymp- 
totically flat spacetime or not, although it is hard to imagine an asymptotically 
flat spacetime filled with elastic material. In any case, it would be nice to be able 




and 
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to tell, by looking at a given constitutive relation, whether the solution under that 
constitutive assumption represents an elastic ball of finite radius or infinite radius. 
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Appendix 



In §2, we have assumed that the material metric is flat. Since there is no 
reason to exclude non flat material space, we will consider the case with non-flat 
material metric. It turns out that the arguments used to prove the local existence 
and global properties in §3 and §4 applies to the case of non flat material space as 
well.. 

Because of the spherical symmetry, the most general form allowed for the 
material metric has the form 

dy 2 + f 2 (y) [d6 2 +sm 2 e d0 2 } 

where / is a non negative smooth function such that f(y) = y ■ g(y 2 ) for some 
smooth function g with g(0) = 1. We will impose an extra assumption that /' > 
and g is a bounded function on R. With this material metric, the strain tensor is 

S = -\og^= (e x ^dr 2 ) - log^M ( r W) _ i og M ( r W0# 2 ) 
v / e A(r) V / r r 

Note that this expression can be obtained simply by replacing in (6) y(r) by 
f(y(r)). As in §2, we proceed to choose the invariants u and v. This time we 
choose 

u:=\og\— L=\ (18) 




a + u = -31og (49) 



Then we have the same stress-strain relationship (10) as in the case of flat mate- 
rial metric. Consequently the expressions in the the Einstein equations and the 
Tolman-Oppenheimer-Volkoff equation (12)- (24) are unchanged. The only devia- 
tion from the flat material case is that auxiliary equation (for v ') is different. 

2xv' = 3 - 3Vbe u ■ f(y) (50) 

Here we can consider y as an expression written in terms of x and v since / is 
invert ible near y = 0. 

We will show that the statement of Theorem 2 in §3 about the local existence 
of solutions is also true for the system (26), (27), (50). First, with the substitutions 
u = xu\ and v = v + xv\ , we will show 

x~ l (l - Vbe u ■ f'(yj) =-w- Ul + 0(x). (51) 
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Once this has been shown, the local existence result for the system (26), (27), (50) 
follows immediately by following the proof of Lemma 1. (Recall that in the proof 
of Lemma 1 the auxiliary equation was used only for the verification of (35), 
which in the current case corresponds to (51).) The proof of (51) follows. From 
(49), y = f~ l (re~§) = 0(r) as r — > 0. So, by the assumption on /, we have 
f'( y ) = g (yi) + 2y 2 g' \y 2 ) = 1 + 0(y 2 ) = 1 + 0(r 2 ) = 1 + 0{x). So, the left hand 
side of (51) can be treated as if there were no f'{y) in the second term, which 
is the case in the proof of Lemma 1. Thus (51) has been verified and the local 
existence follows. 

To obtain the global result stated in Theorem 3 in §4, we only need changes 
in the arguments leading to the existence of upper and lower bounds of v. To show 
that v is bounded from above, in §4, we have used the fact that the second term on 
the right hand side of (50) is positive. The same argument applies in the current 
case since we have assumed /' > 0. To show that v is bounded from below, it 
suffices to show that the second term on the right hand side of the equation (50) 
is bounded from above. The upper bound of u and b := e x can be obtained 
independently of the equation (50) (See the proof of Lemma 2 and the Step II 
which follows this lemma.) So, it remains to show that f'(y) is bounded from 
above. To show that f'(y) is bounded, it suffices to show that y is bounded since 
/' is a continuous function. This will be done by using the equation (48) . Since we 
know that e A and u are bounded from above, it follows from (48) that y' < Kf(y) 
for some constant K. But, we have assumed that f(y) = y g(y 2 ) < y ■ C for some 
constant C. So, y' < Ky, and we have the bound for y. 
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